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ABSTRACT 

The ultra-luminous X-ray (ULX) source ESO 243-49 HLX-1, which reaches a maximum luminosity 
of 10 42 erg s _1 (0.2-10 keV), currently provides the strongest evidence for the existence of intermediate 
mass black holes. To study the spectral variability of the source, we conduct an ongoing monitoring 
campaign with the Swift X-ray Telescope, which now spans more than two years. We found that 
HLX-1 showed two fast rise and exponential decay (FRED) type outbursts in the Swift XRT light- 
curve with increases in the count rate of a factor ^40 separated by 375±13 days. We obtained 
new XMM-Newton and Chandra dedicated pointings that were triggered at the lowest and highest 
luminosities, respectively. From spectral fitting, the unabsorbed luminosities rang ed from 1.9 x 10 40 to 
1.25 x 10 42 erg s _1 . We confirm here the detection of spectral state transitions from HLX-1 reminiscent 
of Galactic black hole binaries: at high luminosities, the X-ray spectrum showed a thermal state 
dominated by a disk component with temperatures of 0.26 keV at most, and at low luminosities the 
spectrum is dominated by a hard power law with a photon index in the range 1.4-2.1, consistent 
with a hard state. The source was also observed in a state consistent with the steep power law state, 
with a photon index of ~ 3.5. In the thermal state, the luminosity of the disk component appears to 
scale with the fourth power of the inner disk temperature which supports the presence of an optically 
thick, geometrically thin accretion disk. The low fractional variability (rms of 9 ±9%) in this state also 
suggests the presence of a dominant disk. The spectral changes and long-term variability of the source 
cannot be explained by variations of the beaming angle and are not consistent with the source being 
in a super-Eddington accretion state as is proposed for most ULX sources with lower luminosities. 
All this indicates that HLX-1 is an unusual ULX as it is similar to Galactic black hole binaries, which 
have non-beamed and sub-Eddington emission, but with luminosities 3 orders of magnitude higher. 
In this picture, a lower limit on the mass of the black hole of >9000 M can be derived, and the 
relatively low disk temperature in the thermal state also suggests the presence of an intermediate 
mass black hole of a few 10 3 M Q . 

Subject headings: X-rays: individual (ESO 243-49 HLX-1) - X-rays: binaries - black holes - accretion, 
accretion disks 



1. INTRODUCTION 

At present, only two families of black holes 
have convincing observational evidence: stellar-mass 
black holes detected in som e X-ray binaries (e.g. 
iRemillard fc McClintockl 120061 ) and super-massive black 
holes (SMBHs, 10 6 " 9 M w ) which are ubiquitous i n 
the center of galaxies ()Kormendv fc Richstond 119951 ). 
sometimes revealing themselves as active galactic nu- 
clei (AGN). Remarkably, the existence of intermediate 
mass black holes (IMBHs, 10 2-5 Mq) remains to be 
proven. Such objects may have been ejected follow- 
ing an interaction with the ce ntral SMBH and may re- 
main in the halo of galaxies (jMicic etalll20Ti . They 
may be found in globular clusters, and the strongest 
case is the massive cluster Gl in M31, where dynami- 
cal, X-ray, and radio studies are consistent with a mas s 
of 2 x 10 4 Mq (IGebhardt et all 12002 IKong et alJiMoh . 
However, radio observations of Galactic globular clus- 
ters currently provide only upper limits on the mass (e.g. 



IMaccarone fc Servihatl WM ILu fc Kondl2011h . IMBHs 
may also be the nuclei of satellite galaxies capt ured dur- 
ing hierarchical merging (jKing fc Dehner3l2005l ). Finally, 
they may be the e ngine of some ult ra-luminous X-ray 
sources (ULXs, e.g. IMiller et alJl200l . 

A ULX is a non-nuclear extragalactic X-ray source that 
has an X-ray luminosity that exceeds the Eddington lu- 
minosity — 1.3 x 10 38 (Mbr/Mq) erg s _1 — for a stellar- 
mas s black hole (~ 10 Mq), supposing isotropic emission 
(see lRobertsll2007l for a review). Different solutions have 
been proposed for the high luminosity problem of ULXs: 
(1) they are stellar-mass black holes emi tting up to a fac- 
tor of 10 above their Eddington limit (|Begelmanl 120021 ) 
in a new "ultra-lu minous" state with sup er-Eddington 
accretion rates (e.g. iGladstone et al.ll2009l ). (2) they are 
stellar-mass bl ack holes with geom etric (e.g. lKindl2009T ) 
or relativistic (jKording et al.ll2002f) beam ing, or (3) the 
objects are IMBHs (e.g. IMiller et alJl2004D . It could also 
be a combination of these possibilities. For the most 
luminous ULXs above 10 41 erg s _1 , the mass argument 
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seems to be the main one that can convincingly explain 
the extreme luminosities. 

One of the keys to further constraining the nature 
of the black hole in such sources is the study of 
their spectral and timing variability. Indeed, Galac- 
tic black hole binaries (GBHBs) containing stellar- 
mass black holes have been commonly observed to un- 
dergo transitions between different spectral states (e.g. 
i Tananbaum et al.lll972t iKubota k Makishimall2004 see 
Remillard k McClintockll2006t iDone et al.ll2007t iBellonil 
2010 for reviews). In the thermal state, the emission is 
dominated by an optica lly thick, geometrically th in ac- 
cretion disk component (iShakura k Sunva cv 1973) , with 
temp eratures of ~1 keV (e.g. iRemillard k McClintockl 
2006). For such a disk extending down to the in- 
nermost stable circular orbit (ISCO), the luminosity 
Ldisk, inner temperature Tj n and mass of the black 
hole Mbh are theoretically related: idisk °c for a 

given M B h, and T in oc M b "h /4 (|Shakura k SunvaevlU97l 
iMakishima et all [2000) . In the hard state, the accre- 
tion disk appears to be fainter and cooler, and may be 
truncated at a large radius. The physical condition of 
material within this radius remains uncertain. Investi- 
gations of GBHBs in the hard state suggest that both 
synchrotron and Compton components contribute to the 
broadband, power law like spec trum with typical photon 
indic es of 1.4 < T < 2.1 (e.g. IRemillard k McClintockl 
2006). Quasi-steady radio jets are commonly observed 
during this state, and clear correlations bet ween the radio 
and X -ray intensities have been reported ([Fender et al.1 
I2004T) . Some GBHBs also showed a steep pow er law state 
(r > 2.4, e.g. IRemillard k McClintockl 120061) . 

Consideration of simple properties such as X-ray spec- 
tral hardness and fractional variability have led to a 
better u nderstanding of the physical processes in ac- 
tion (e.g. [Bclloni 2010). In a hardness-intensity diagram 
(HID), the source shows a hyster esis curve between the 
thermal state and the har d state ([Miyamoto et al.lll995t 
iMaccarone k Coppill200l . The root mean square (rms) 
variability amplitude is generally low in the thermal state 
(<10%), and can reach higher values in the hard state, 
with typical rms of ~30%. As seen in a hardness-rms 
diagram (HRD), the rms is anti-corr elated with fl ux and 
positively correlated with hardness (Bclloni 2010). 

The lack of a canonical thermal state, with a dominant 
disk component , seems to be a common feature of ULXs 
(IFeng k Kaareti [20051 : iWinter et all [2001 iSoria et al.1 
2009). Two ULXs showed such a thermally dominated 
spectrum: M82 X41.4+6 0=X-1 (kT ln = 1.1-1 5 keV and 
L x = 2-8 x 10 40 erg s' 1 . IFeng k Kaaredl2010h . and M82 
X37.8+5 4 (0.9 -1.6 keV and L X max = 4.4 x 10 39 erg s" 1 , 
Uin et alJl2010h . Those sources could harbor a fast spin- 
ning black hole of 200-800 M and <100 M , respec- 
tively. ULXs with fcTjn ~ 0.1 keV soft disk spectra 
may be good candidates for IMBHs in the thermal state. 
Some ULXs have shown soft excesse s that could be fit - 
ted by such a cool disk emission (e.g. lMiller et"aLl l2003). 
However, this component is not dominant, thus the 
sources are not in the canonical thermal state. A few in- 
triguing ultraluminou s supersoft sources show tempera- 
tures of ~0.1 keV (e.g.lFabbiano et al.ll2003t iMukai et al.1 
120031: iKong et al J 120041 Uin et all 120 111) , but they could 
be explained by matter outflow from super-Eddington 



accretion (King k Poundsl[200l . 

ULXs do not generally show quali tative spectral 
chang es when the luminos ity varies (IFeng k Kaareti 
[20051 iKaaret k Fend [20091 : IGrise et all[2010j) . PossT 
ble state transitions have been reported for some ULXs 
but they are subtle and not c learly similar to GB- 
HBs: e.g. Holmberg IX X-l (|La Parola et all 120011: 
Vierdavanti et all [2010) or M82 X37.8+54 ((Jin_et^l| 



2010). 



ULXs are generally highly variable sources show- 
ing variations of fluxes by a factor of ~5-15 on 
time scales of hours to yea rs (|La Parola et all 120011 : 
IGrise et all 120101: |Kong|l201lD . On shorter time scales, 
six ULXs showed intrinsic variability with power spec- 
tra in the form of e ither a power-la w or broken-power- 
law-like continuum (|Heil et all 120091) and in some cases 
quasi-periodic oscillations (QP Os, iCasella et al.l 120081 : 
IStrohmayer k Mushotzkvll2009h . In some ULXs the vari- 
ability is signific antly suppressed compared to bright GB- 
HBs and AGN (IHeil et alJl2009h . 

To date, one of the best IMBH candidates is the most 
luminous ULX ESO 243-49 HLX-1 (HLX -1 for short) 
which reaches a luminosity of 10 42 erg s _1 (jFarrell et all 
2009). The distance, and thus the high X-ray luminosity 
of HLX-1, have been firmly confirmed through the de- 
tection of an Ha li ne at a redshift consi stent with that 
of its host galaxy (jWiersema et~^[2010l ). The Ha line 
has been found in the spectrum of the optical counter- 
part which falls inside the Chandra error circle of the 
source ( 0"3 at a 95% con fidence level iWebb etaLl [20101: 
see also ISoria et al.lfeOlOl ). Observations of HLX-1 with 
the Swift X-ray Telescope (XRT) in August 2009 showed 
a hardening of the source at low luminosities which was 
interpreted as the first evidence f or a change to the hard 
spectral state ([Godet et alll2009T ). 

Our aim in this paper is to measure the spectral hard- 
ness and the variability of HLX-1 at different luminosity 
states, using all the X-ray data available from current ob- 
servatories!]] in order to test similarities with GBHBs or 
other ULXs. For this purpose, we obtained new XMM- 
Newton and Chandra dedicated pointings that were trig- 
gered at the lowest and highest luminosities, respectively. 
The X-ray datasets are described in Section [2l Spectral 
analysis of firstly the high luminosity states (Section [3]) 
and then the low luminosity state (Section [4j follows. A 
global analysis is given in Section [5] and the final Sec- 
tion [5] discusses the properties of HLX-1 compared to 
GBHBs and other ULXs. 

2. DATA AND LONG TERM VARIABILITY 

2.1. Swift XRT data 

HLX-1 has been regularly monitored by the Swift XRT 
since October 2008. Since August 2009, observations oc- 
curred on average every ~l-2 weeks. All the data were 
processed with the Swift XRT pipelin43 version 0.12.4. 
We used the grade 0-12 events, giving slightly higher ef- 
fective area at higher energies than grade events. The 
background extraction region, chosen to be close to the 
source extraction region, is the same for all epochs. No 
XMM-Newton sources are present inside the background 



1 Except two Chandra HRC-I observations (Webb ct al. 2010) 
as this instrument does not provide spectral information. 

2 http:/ /heasarc. gsfc.nasa.gov/docs/swift/analysis/ 
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Figure 1. 

dates of the 



Swift XRT light curve of HLX-1. The epochs that were grouped to generate spectra are shaded with labels SI to S10. The 
XMM-Newton and Chandra observations are represented by vertical lines. 



Table 1 

XMM-Newton observations of HLX-1 



Obs. name 


ObsID 


Date 


MJD 


Exp. time [ks] 


Instrument 


Mode 


Source region 


Background region 


XMM1 


0560180901 


2004-09-23 


53271 


21 


MOS1&2, pn 


full-frame 


22.5" 


50.0"-63.4" 


XMM2 


0560180901 


2008-11-28 


54798 


50 


MOS1&2 


full-frame 


27" 


60" -76" 












pn 


small window 


27" 


3x27" 


XMM3 


0655510201 


2010-05-14 


55330 


100 


MOS1&2 


full-frame 


14" 


40"-47" 












pn 


small window 


14" 


3x14" 



extraction region. 

We generated a light curve from all the Swift point- 
ings with a binning of a minimum of 50 counts per bin 
for the en ergy band 0.3-10 keV (F igure [TJ using the web 
interface (|Evans et al.ll2007L l2009h . There is a recurrence 
in the light curve of 375±13 days between the two in- 
creases that may be periodic, but this will have to be 
confirmed by further observations. If the behavior is in- 
deed periodic, the next maximum would occur in Au- 
gust/September 2011. The cause o f this variabi l ity ha s 
been investigated in more details bv lLasota et al.1 (|2011h . 

We combined the data into different epochs labelled 
SI to S10 in Figure [TJ and extracted a spectrum for 
each epoch. A 20-pixels (47'.'2) radius circle was used 
to extract the source and the background spectra us- 
ing XSELECT v2.4a. The ancillary files were created 
with XRTMKARF vO.5.6 and exposure maps gener- 
ated with XRTEXPOMAP vO.2.7. The response file 
swxpc0tol2s6_20070901v011.rmf was used in the spec- 
tral fitting process. We used a binning of a minimum of 
20 counts per bin for each spectrum except the S2 and 
S6 sp ectra in order to use the x 2 statistic within Xspec 
12.6 (|Arnaudlll996ri. For the S 2 and S6 spectra, we used 



the Cash statistic ()Cashlll979l ) due to the lower number 
of counts in these spectra. 



2.2. XMM-Newton EPIC data 

The field of ESO 243-49 and HLX-1 was observed three 
times with XMM-Newton (TableUJ. HLX-1 was detected 
serendipitously at an off-axis angle of 9.'29 in XMM1, and 
XMM2 was obtained to sh ow the spectral variability of 
HLX-1 (|Farrell et al.ll2009r i. Following a transition to a 
very low flux state, we triggered the XMM3 observation. 

The data were processed using the XMM-Newton Sci- 
ence Analysis Systen0 (SAS) vlO.O software with the 
most recent calibration files as of 2010 December 13. The 
observation data files (ODF) were reduced using the ep- 
proc and emproc tasks to produce event lists. Single 
event light curves with energies exceeding 10 keV were 
generated for each camera in each observation in order 
to identify periods of high background related to soft 
proton flares. In the XMM1 and XMM2 observations 
the background levels were low in each camera with no 
flaring events, and so no good time interval (GTI) filter- 
ing was applied. However, during the XMM3 observa- 
tion significant flaring events were present at the start 
and end of the exposure, with an additional small short- 
duration flare occurring ~16 ks into the observation. For 
the XMM3 data we therefore generated GTI files for the 

3 http://xmm.esa.int/sas/ 
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Figure 2. Hardness intensity diagram for HLX-1 based on XMM- 
Newton count rates (counts s — 1 ). Hardness ratios are HR = 
(B2-B1) / (B1+B2) with Bl: 0.3-1 keV and B2: 1-10 keV. The 
discrepency between MOS1 and 2 in XMM-1 was already identified 
in Farrell et al. (2009). 



spectral extraction using the SAS task tabgtigen and cut- 
off count rates of 0.3 counts s _1 and 0.1 counts s _1 in the 
flare background light curves for the MOS and pn cam- 
eras respectively. This filtering resulted in net exposure 
times of ^98 ks and ^99 ks (for pn and MOS respec- 
tively). The same event filtering criteria as used for the 
production of pipeline products for the 2XMM catalogue 
were used to produce light curves and spectra for each 
camera in each observation (e.g. single to double events 
for the pn and sing le to quadruple events for the MOS, 
iWatson et alll^OOl . 

Source and background spectra were extracted for each 
camera for each observation, with response and ancillary 
files generated in turn using the tasks rmfgen and arfgen. 
The source extraction region radii (reported in Table [lj 
were chosen so as to optimize the signal-to-noise based on 
the detected count rates and off- axis position of HLX-1 
in each of the observations (V. Braito, private communi- 
cation) . The background extraction regions were chosen 
to be three times the area of the source extraction re- 
gions, so as to provide a robust estimate of the average 
background level at the position of HLX-1. The spec- 
tra were grouped to have at least 20 counts per bin to 
provide sufficient statistics for spectral analyses using x 2 
statistic. 

Source light curves were extracted using events in the 
energy range of 0.3-2 keV for the pn camera of all obser- 
vations. They were binned at the frame time, which is 
73.4 ms for XMM1 and is 5.9 ms for XMM2 and XMM3. 

Figure [2] shows the count rates for each camera as a 
function of a hardness ratio. We see a discrepancy be- 
tween the XMM1 MOS1 and MOS2 hardness ratios, with 
the MOS2 data showing a deficit in counts in the 0.3- 
1 keV energy band compared to both the MOS1 and 



pn dat£0- We thus excluded energies below 0.5 keV 
for MOS2. Despite this inconsistency between the MOS 
hardness ratios in the XMM1 data, a clear trend is ob- 
served indicating significant spectral hardening at lower 
luminosities, which is most clearly shown by the pn data 
(Figure 2, top panel). 

2.3. Chandra ACIS-S data 

We obtained a 10 ks observation with the AXAF CCD 
Imaging Spectrometer (ACIS) onboard Chandra on 2010 
September 6 (ObsID 13122, MJD 55445) in the director's 
discretionary time prog ram after the sourc e showed a re- 
brightening in X-rays (|Godet et al.ll2010D . In order to 
limit pile-up, this observation was performed with the 
ACIS-S3 chip only and using a 1/4 chip subarray with a 
frame time of 0.8 s, leading to 5% dead time (0.041 s). 
The ACIS-S3 chip was chosen for its higher sensitivity to 
softer X-ray photont0. 

The data were processed u sing CIAC04.2and CALDB 
4.3.0 dFruscione et all 120061) . We recreated the evt2 file 
following the corresponding thread^ to apply the latest 
calibration files. We extracted the source spectrum using 
psextract with an extraction circle of radius 2'.'2 which 
includes 90% of the energy at 8 keV. The background 
spectrum was extracted in an annulus surrounding the 
source with inner and outer radii of 15" and 50" respec- 
tively (thus excluding the host galaxy). We then ran the 
tasks mkacisrmf and mkarf to obtain the instrument re- 
sponse files as recommended for ACIS-S data. The spec- 
trum was binned with a minimum of 15 counts, and we 
discarded bins below 0.3 keV where the response of the 
instrument is not calibrated. 

3. HLX-1 AT HIGH LUMINOSITIES 

3.1. Re-analysis of XMM1 and XMM2 spectra 

Spectral fits of the XMM2 data with th e relativistic 
model BHSPEC have been performed by iDavis et alJ 
(|2011h . and a more detailed analysis of all spectra with 
physically motivated accretion disk models used for GB- 
HBs and ULXs is presented by Godet et al. (2011, in 
preparation). We used here simple models to fit the X- 
ray spectra, e.g. an absorption on the line of sight, a 
power law model (po w), and a multi-temper ature black 
body model (diskbb, iMakishim a et al.lll986f). We used 
the wabs absorption model ([Morrison fc McCammonl 
1983). The Galactic absorption along t he line of sight 
of ES Q 243-49 is 1.8 x 10 20 atom cm" 2 (jKalberla et all 
l2005h . 

For XMM1, we fitted simultaneously the MOS1, MOS2 
and pn spectra using Xspec and ignored channels with 
energy lower than 0.3 keV as well as channels tagged as 
bad. The spectra are well represented by a simple ab- 
sorbed steep power law of index 3.5^^2 g ivm g a reduced 
X 2 of 0.93. The addition of a diskbb component improved 
the x 2 by 3.4. To test the significance of this compo- 
nent we used the Bayesia n posterior pred i ctive proba- 
bility values method (e.g. iProtassov et al.l I200 2). as it 

4 This discrepancy had been noted earlier by Farrell et al. (2009), 
leading them to exclude energies below 0.5 keV for MOS2. 

5 http://cxc.harvard.edu/proposer/POG, Proposers' Observa- 
tory Guide 

6 http://cxc.harvard.edu/ciao/ 

7 http:/ /cxc. harvard.edu/ciao/threads/createL2 
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Table 2 

Parameters of best fits for HLX-1 X-ray spectra 



Obs 
(1) 


-Wh20 
(2) 


kT in 
(3) 


r 

(4) 


X 2 /dof 
(5) 


-^[0.2-10 keV] 

(6) 


unabs. Lr . 2 _ 10 keV] 
(7) 


unabs. L disk 
(8) 


Steep power law state 
XMM1 10 ± 4 
XMM1 5 ± 5 


27+ 12 
U - Z '-0.05 


o c+0.3 
J.O_ Q 2 

o 9 +0.3 


59.7/64 
56.3/62 


3.2 ±0.5 
3.1 ± 1.2 


9.9 ±0.5 
6.9 ±2.6 


1 6 +0 ' 7 

l.u_ 1 2 


Thermal state 
















Chandra^ 

XMM2 
S1+S4+S8 

S3 

S7 


2+ 4 
3±1 

3 

3 

3 


0.22 ± 0.02 
0.18 ± 0.01 
0.18 ± 0.02 
0.26 ± 0.03 
0.23 ± 0.03 


2.1 ± 0.2 
9 Q+0.7 
z -°-0.5 


57.0/48 
395.6/346 
48/45 
14.1/21 
12.3/14 


6.7 ±0.3 
3.5 ± 0.2 

8.1 ±0.5 
7.7 ±0.5 


9.0 ±0.3 
6.4 ± 0.2 

12.2 ± 0.5 
12.5 ± 0.5 


13.0 ±0.3 
6.6 ± 0.3 

15.2 ±0.5 
16.5 ±0.5 


Hard state 
XMM3 
XMM3 
XMM3 b 
XMM3 b 
S2+S6 


6±S 
62if 
3 
3 
3 


0.09 ± 0.02 
0.07 ±0.04 


9 c+0.6 

9 4 +0.5 
z '^-0.8 

2.0 ±0.3 
1.6 ±0.4 

2.1 ±0.4 


85.2/60 
69.3/58 
63.6/59 
58.0/57 
8.4/8 


17+ - 01 

-0.09 

17+ - 02 
U - i '-0.05 

0.16 ±0.03 
0.20 ±0.04 
0.6 ±0.2 


31+ - 02 

~21 
0.19 ±0.03 
0.28 ± 0.04 

0.8 ±0.3 


~56 
0.4 ±0.1 



Note. — Columns: (1) Observations; (2) Absorption on the line of sight, when no error, it is fixed to the best 
estimate of 3 ± 1 X 10 20 atom cm" 2 for XMM2; (3) Temperature of the diskbb model in keV and 90% error; (4) 
Photon index of the power law component and 90% error; (5) \ 2 °f the fit and degrees of freedom; (6) Absorbed flux 
in the 0.2—10 keV energy range for the model, with la error in 10 — 13 erg s — 1 cm -2 ; (7) Unabsorbed luminosity in the 
0.2-10 keV energy range assuming a distance of 95 Mpc, with lcr error in 10 41 erg s — x ; (8) Unabsorbed bolometric 
luminosity of the disk component assuming a distance of 95 Mpc, with lcr error in 10 41 erg s — 1 . 
a The Chandra spectrum was fitted with a pile-up model with a frame time of 0.8 s, alpha=0.99 and psffrac=0.95, 
which was then removed to estimate fluxes. 

b We added a mekal component at a fixed redshift of z = 0.0224 and with solar abundances which significantly 
improved the fit. We found a temperature kT = 0.43 ± 0.09 keV. This component was removed before estimating 
the flux and luminosity of HLX-1, assuming it is due to the host galaxy. 



was already done by lFarrell et al.l (|2009| ). We generated 
3000 simulated set of spectra with the Xspec fakeit com- 
mand using the absorbed power law model and the re- 
sponse files corresponding to the XMM1 spectra. We 
then computed the distribution of the x 2 improvement 
of the fit when adding a diskbb component. The x 2 i m ~ 
provement of 3.4 for XMM1 translates to a significance 
of 73% for the diskbb comp onent, insufficient to claim it 
is real (|Farrell et al.l 12003 derived a significance of 70% 
with a similar method). 

For the XMM2 spectra, we included channels with en- 
ergy higher than 0.2 keV due to the higher signal to noise 
at low energies. The spectra are poorly fitted with a sim- 
ple absorbed power law mode (x 2 of 577.7 and 348 de- 
grees of freedom) . When adding a diskbb component the 
reduced x 2 significantly dropped to an acceptable value 
of 1.14, indicating that this model better represents the 
spectra. 

For XMM1 a nd XMM2, we thus obtained best fits in 
agreement with lFarrell et all ()2009f ). We show the corre- 
sponding folded spectra in Figure |U The calibration was 
improve d in the low energy range with respect to that 
used by iFarrell et al.l (|2009T ) which lead to smaller error 
bars in this work. We thus present new estimates of the 
parameters in Table [2] 

3.2. Chandra spectrum and pile-up 

The spectrum was fitted in Xspec with an absorbed 
diskbb model which gave a reduced x 2 of 1.43 for 49 
degrees of freedom. A clear hard excess appeared in 
the residuals th at could correspond to the effect of pile- 
up ()Davisl l200lT) . We indeed estimated a possible level 



of pile-up of 5% with the Chandra Proposal Planning 
Toolkit^ ba sed on PIMM S v4.2. We thus added a pileup 
component (Davis 2001) to the model and obtained a sat- 
isfactory fit with a single absorbed diskbb model with a 
reduced x 2 of 1.18 (see Figuregland Tabled]). We let the 
parameter alpha free and found a value of l.O+p'g. This 
parameter is related to the probability of events being 
retained as a good grade after filtering and we would ex- 
pect a value between 0.5 and 0.7. We note that the value 
of the diskbb temperature is not significantly affected by 
this parameter. Large contributions to the x 2 seem to 
come from features around 0.6 keV (emission) and 1.1 
and 1.3 keV (absorption). However, given the low num- 
ber of degrees of freedom, there is a large uncertainty on 
the expectation of the reduced x 2 b eing 1 (a of 0.2 on the 
estimate of the reduced x 2 , see e.g. Andrae 2010), so we 
cannot claim that these features are real. We therefore 
consider the continuum well fitted by the diskbb model. 
We found that an additional power law component with 
a fixed photon index of 2 would contribute to at most 
10% of the 0.2-10 keV absorbed flux, indicating that the 
disk component is dominant in the spectrum. 

The Chandra observation was simulated to better char- 
acterize the emission of HLX-1. We used the ray tracing 
tool ChaRT dedicated to Chandra (jCarter et al.l 120031) 
and ge nerated an image with MARX for the ACIS-S de- 
tector (jWise et al.lll997D . The best fit spectral model of 
the source, obtained as explained above, was used as an 
input in the energy range 0.3-8 keV. A simulated image 
is presented in Figure [5j 



http: / /cxc. harvard, edu /toolkit/ pimms . j sp 
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Figure 3. XMM-Newton folded spectra of HLX-1 and best fit 
models. Data points from the three cameras are reported with 
their error bars and the model is indicated as a dashed line. Dotted 
lines show components of the model. 
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pileup.abs(diskbb) 
xVdof = 57.0/48 




Energy (keV) 

Figure 4. Chandra ACIS-S folded spectrum of HLX-1. Best fit 
(top) and residuals (bottom). See Table [2] 

In the band 0.3-8 keV, the results from the simulation 
gave 1294 ± 36 counts. After running the pileup tool in 
MARX, we obtained 1182 ± 34 expected counts (9 ± 5% 
loss). We detect 1130±34 counts (14±5% loss), which is 
consistent within the errors with the MARX simulation 
with pile-up effects. The PSF in the observation shows 
a deficit of counts in the central pixels compared to the 
simulation without pile-up. It is however consistent with 



HLX-1 



Chandra ACIS-S 




ChaRT + MARX simulation 



Figure 5. Chandra ACIS-S images of HLX-1 in the 0.3-8 keV en- 
ergy band. Top: data. Bottom: result of simulation using ChaRT 
and MARX with pile-up included. Galaxy contours are overlaid 
(at H-band Vega mag arcsec -1 of 21, 18.5 and 16). A dashed 
black 5" radius circle is shown around the bulge of the galaxy and 
similar grey circles are placed around the galaxy for comparison. 
The position of HLX-1 is indicated with its 0"3 error at 95%. 

the image of the simulation with pile-up applied. We are 
therefore confident that the use of the pile-up model in 
our spectral fitting is justified. 

3.3. Emission from the galaxy 

No significant point source is detected inside the galaxy 
contours except HLX-1. We looked for an excess of 
counts in an extended 5" radius region encircling the 
bulge of the galaxy. In the 0.3-8 keV energy band, we 
expect 5 ± 2 background counts in such a region and we 
may have in addition 5 ± 2 counts from HLX-1 due to 
the scattering of its emission on the detector (Figure [SJ 
lower panel). Finally, we expect a maximum of 10 ± 3 
counts, while 19 counts were detected. We repeated the 
same estimation for different positions on the detector 
and found no other excess of counts (Figure 03 upper 
panel). We found a similar excess in the 0.5-6 keV energy 
range where the signal-to-noise is the best for ACIS-S: 15 
counts detected, background of 3 ± 2 counts and HLX-1 
emission of 3 ± 2 counts. 

The median energy from those photons is 0.92 ± 
0.10 keV suggesting a soft emission. The emission from 
the bulge of spiral galaxies have previously been mod- 
eled by a mekal model with tem peratures of 0.3 to 0.6 
keV (e.g. Humphr ey et al.l I2004J ). Assuming a mekal 
model (at the redshift of the galaxy, z = 0.0224, with 
solar metallicity and fixed absorption to the Galactic 
value of 1.8 x 10 20 atom cm -2 ), we find a temperature 
kT = 0.30 ± 0.15 and a flux of - 5 x 10~ 15 erg s" 1 cm~ 2 
in the 0.2-10 keV range. At the distance of the galaxy, 
this converts to an unabsorbed X-ray luminosity of ~ 
6 x 10 39 erg s _1 . This level of luminosity is consistent 
with the integrated X-ray l uminosity of spir al galaxies in 
general (10 38 -10 42 erg s" 1 . lFabbianolll989f) . This emis- 
sion might contaminate XMM-Newton and Swift HLX-1 
data at low luminosities given their lower angular reso- 
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Figure 6. XMM3 spectrum residuals for different models reported 
in Table [2] 

lution. 

4. HLX-1 AT LOW LUMINOSITIES 
4.1. XMM- Newton spectrum XMM3 

We fitted simultaneously the MOS1, MOS2 and pn 
spectra using Xspec and ignored channels with energy 
lower than 0.3 keV as well as channels tagged as bad. We 
used an absorbed power law and obtained a reduced x 2 
of 1.42 with 60 degrees of freedom (Tabled]). The high- 
est contributions to the x 2 are associated with a possible 
soft excess around 0.7 keV and a hard excess at high en- 
ergies (Figure HI first panel). We thus tested two possible 
additive components to the model to account for those 
features: (1) a diskbb component to test for the presence 
of a thermal disk, and (2) a mekal component which 
would correspond to a possible contamination from the 
galaxy bulge to the spectrum (see Section |3~3"|) . 

When adding a diskbb component to the single ab- 
sorbed power law model, the reduced x 2 is lowered 
to 1.20. However, a marked soft excess appears in 
the residuals and the hard excess is still present (Fig- 
ure [SJ second panel). The best fit is found for an 
Nn ^6.5±1.5xl0 21 atom cm -2 , an order of magnitude 
higher than the XMM2 value and a lower temperature 
(kTi n ~ 0.08 keV) than for the Chandra data in the high 
state. The bolometric unabsorbed luminosity of the disk 
component is thus higher than in the high luminosity 
state (Section [3]). 

The addition of an abs(mekal) component — at a fixed 
redshift of z = 0.0224, with solar abundances and Galac- 
tic absorption as it could not be constrained — to the 
single absorbed power law improved the reduced \ 2 to 
an acceptable value of 1.08 and resulted in a lower pho- 
ton index for the power law (2.0 ± 0.3). The absorption 
of HLX-1 is not well constrained, so we fixed it to the 
XMM2 value which is the most precise estimate we ob- 
tained (see Table [2]). We note that if the absorption 
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Figure 7. XMM-Newton MOS images of HLX-1 in the energy 
band 0.5-6 keV. The pixel size is 3'.'2. The H-band galaxy contours 
are shown as in Figure[5] The XMM-Newton position error (<3'.'2 
at 95%) is represented by dashed circles and the Chandra error 
circle (0'.'3 at 95%) is indicated by a red arrow. 

is thawed, it tends to zero while the power law index 
further decreases. The residuals are shown in Figure [5] 
(third panel) and the folded spectrum is shown in Fig- 
ured] We found a temperature of kT = 0.43 ± 0.09 keV 
consistent with the expected emission from the galaxy 
(Section 153]) . Moreover, the flux found for this absorbed 
mekal component is 3.3±1.3x 10 -15 erg s -1 cm~ 2 which 
is comparable to the estimate of the possible galaxy con- 
tribution found in the Chandra data (Section 13.31) . We 
are thus confident that some emission from the galaxy 
contaminates the XMM3 spectrum of HLX-1, at a level 
of ~17% in flux at most. 

Finally, we tested the presence of a diskbb component 
in this last model. The resulting reduced x 2 is 1.02, the 
power law is harder (1.6±0.4) and the residuals seem fea- 
tureless for this fit (Figure IB] fourth panel) . We tested 
the significance of the diskbb using the Bayesian poste- 
rior predictive probability values method. In the same 
way as in Section 13.11 we generated 3000 simulated set 
of spectra based on the absorbed power law plus mekal 
model and the XMM3 response files. We then computed 
the distribution of the \ 2 improvement of the fit when 
adding a diskbb component. The x 2 improvement of 5.6 
for XMM3 translates to a significance of 93% for the 
diskbb component. We note that a diskbb component is 
not needed to adequately describe the spectrum, but this 
possibility is intriguing and cannot be ruled out. 

4.2. Astrometry consistency 

We aligned all images on the Chandra image, taking as 
a reference 3 point-like X-ray sources surrounding HLX- 
1. This allowed us to reduce the positional error by 
replacing the XMM-Newton Absolute Measurement Ac- 
curacy of 4" by a residual error of less than 1" after 



§4.7 



the XMM-Newton Users Handbook, 
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alignment. We still have to take into account the Rela- 
tive Pointing Error of 2" and the error on the centroid 
of the source, which yielded an uncertainty lower than 
the size of the pixels (3'.'2) for all the images presented 
in Figure [7J We also note that the relative astrometry 
within all EPIC cameras is accurate to better than 1'.'5. 
All errors are 95% confidence. 

We can expect that ESO 243-49 contributes to the 
X-ray emission of HLX-1 in XMM3, as determined us- 
ing the higher angular resolution Chandra data (Sec- 
tion [33]), since the bulge diffuse emission should be con- 
stant. The source is slightly offset from the Chandra 
position which falls just outside the 95% error circle. It 
is shifted towards the center of the galaxy ESO 243-49, 
indeed suggesting a contamination from the bulge of the 
galaxy. However, the source detection task emldetect did 
not split the source into two sources, indicating that the 
statistic is not sufficient to claim the presence of a con- 
taminating source in the imaging data such as the galaxy 
bulge. The PSF full width half maximum (FWHM) is 5" 
and 6" for MOS and pn, respectively, and the half en- 
ergy width (HEW) is 14" and 15". The distance between 
HLX-1 and the galaxy center is about 8" so the images 
cannot be used as evidence for contamination. 

5. HARDNESS-INTENSITY AND HARDNESS-RMS 
DIAGRAMS 

5.1. Swift XRT spectra 

We completed the study of the Chandra and XMM- 
Newton spectra with the spectral fitting of grouped spec- 
tra from Swift XRT (Figure [TJ. We fitted each Swift 
XRT spectrum using Xspec. For all spectra, given the 
reduced statistical quality, we performed a basic fitting of 
the data using the XMM2 model (i.e. abs(diskbb+pow), 
see Table [2]) multiplied by a constant factor in order 
to estimate the unabsorbed luminosity. For each spec- 
trum we estimated the hardness as the ratio of the count 
rates in the 1-10 and 0.3-1 keV bands. We then used 
WebPIMMS^I to convert those ratios into flux ratios: in 
the low and high states (power law with photon index 
~ 2.1 and ~ 4.0, respectively) we found a conversion fac- 
tor of 1.5 and 0.9, respectively. We also used the follow- 
ing model: abs(Clxdiskbb + C2xpow) where CI and C2 
are two constants. The spectral parameters of the diskbb 
and power law models were fixed at the values found in 
XMM2. The goal is to investigate the evolution in the 
contribution between the soft and hard components. All 
results are reported in Table [3] 

For the S3 and S7 spectra (high luminosity states, see 
Figured]), we used a simple absorbed diskbb model and 
reported the parameters of the fits in Table [5] There is 
no evidence in these spectra for the need of a hard energy 
component. We note that the temperature for S7 is con- 
sistent with that derived from the Chandra observation 
which was performed at a similar time. The normaliza- 
tion is lower for the Chandra spectrum, but the fluxes 
for both spectra remain consistent within the errors. 

The basic spectral fitting suggests that the source is in 
the same state in SI, S4 and S8. We merged the data 
from SI, S4 and S8 that approximately overlap in lumi- 
nosity assuming a ~375 days periodicity of the source. In 
order to get a good fit, it is necessary to add a hard com- 

10 http:/ /heasarc. nasa.gov/Tools/w3pimms. html 
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Figure 8. The PDS of XMM-Newton observations of HLX-1 using 
the pn camera. The black solid constant lines are the average PDS 
above 1 Hz, representing the Poisson level. 

ponent (power law) to the model. Results are reported 
in Table [2j The resulting parameters are consistent with 
those found for XMM2 at a similar level of luminosity. 

Finally, we combined S2 and S6 to increase the signal- 
to-noise of the Swift XRT data in the low luminosity 
state. We obtained a good fit with an absorbed power 
law (r = 2.1 ± 0.4, Tabled]). We note that the flux is 
higher than for the XMM3 data. However, as the S2+S6 
spectrum is from data obtained at different epochs and 
integrated over a long time, the comparison with XMM3 
is not straightforward and differences in the two spectra 
are expected. The S2+S6 spectrum may also be affected 
by contamination from the galaxy bulge (Section 13. 3p . 

5.2. Timing analysis 

We calculated the power density spectra (PDS) for the 
light curves of HLX-1 from the three XMM-Newton ob- 
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Table 3 

Fit results for Swift XRT spectra 



Name 


Start 


End 


Hardness 


Count rate 


Unabs. _Lx 


X / d °t 


C 1 


K I 


X /dot 


(i) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


SI 


54763 


54784 


0.34 ±0.04 


1.20 ±0.06 


6.4 ±0.6 


21.2/17 
30/35* 


0.9 ±0.1 


1.3 ±0.4 


19.4/16 
29.5/35* 


S2 


55048 


55050 


2.4 ± 1.2 


0.15 ± 0.09 


0.5 ±0.3 


<0.09 


0.4 ±0.3 


S3 


55059 


55063 


0.34 ± 0.03 


3.30 ±0.13 


11 ± 1 


14.4/21 






Si 


55106 


55164 


0.27 ± 0.03 


1.86 ±0.09 


9.0 ±0.6 


26.6/19 


1.4 ±0.2 


1.6 ±0.5 


25.5/18 


S5 


55170 


55225 


0.34 ± 0.05 


1.12 ±0.06 


5.8 ±0.6 


15.8/12 


0.7 ±0.1 


1.6 ±0.5 


10.0/11 


S6 


55237 


55421 


1.7 ± 1.0 


0.12 ± 0.03 


0.2 ±0.1 


70/63* 





3+ 01 

U.J_ 2 


65.4/63* 


S7 


55437 


55453 


0.21 ±0.02 


2.80 ±0.10 


11.8 ± 0.3 


20.8/24 
8.6/8 






S8 


55462 


55490 


0.21 ± 0.03 


2.00 ±0.10 


9.6 ±0.6 








S9 


55498 


55543 


0.26 ±0.04 


1.52 ± 0.09 


8.3 ±0.6 


11.2/11 


1.3 ±0.2 


1.5 ±0.6 


11.0/10 


S10 


55548 


55562 


0.32 ±0.04 


0.8 ± 0.3 


5.9 ±0.6 


12.6/14 


0.84 ±0.14 


1.2 ±0.4 


10.8/13 



Note. — Columns: (1,2,3) Name, start and end day (MJD) of observation as reported in Figure[T] (4) Hardness 
is given as the ratio of fluxes F[l - 10 keV]/F[0.3 - 1 keV]]; (5) Count rate (10 -2 count s _1 ) in the band 0.3- 
10 keV; (6) Unabsorbed luminosity in 10 41 erg s — 1 after fitting the XMM2 model multiplied by a constant; (7) 
X 2 and degrees of freedom (dof) for the fit; (8,9,10) Result of the fit with the two XMM2 components multiplied 
by factors CI and C2, and goodness of the fit. 

* Cash statistic was used and the number of PHA bins is indicated instead of dof. 



servations. The light curves were split into four segments 
of equal lengths, and the PDS were calculated for each 
segment. For each light curve, all four PDS were merged 
and averaged by binning in frequency using a logarithmic 
factor 1.1, under the condition that each bin contains at 
least 20 individual PDS measurements. The errors were 
calculated from the sample standard deviation of PDS 
measurements in each bin. 

The PDS of HLX-1 are shown in Figure H The black 
solid lines denote the average PDS above 1 Hz. Those 
values deviate from the expected Poisson level by less 
than 0.4%. We see that all PDS are flat, showing no 
significant intrinsic source variability, and no QPO. The 
fit results using a constant Cp are given in Table |U To 
obtain a constraint o n the variab i lity, w e fo llowed the 
proced ure adopted bv lGoad et al.l (|2006l ) and iHeil et al.l 
(2009). In this procedure, the PDS are fitted with a 
broken power-law (BPL) model or a Lorentzian model 
plus a constant. 

For the BPL, the indices below and above the break 
frequency (/b) were assumed to be —1 and —2, respec- 
tively. The value Ci// corresponds to the power times 
the frequency below /b from this model. The upper limits 
of C\/f for two values of /b (10 -3 and 1 Hz) are shown 
in Table |4j We found that these limits are relatively 
independent of the assumed /b value, at least within 
the 10 _3 -1 Hz range. Similarly, it is found to be fairly 
const ant across AGN and GBHBs in the range 0.005- 
0.03 (|Papadakisll2t)0l . The upper limits of C 1/f for 
XMM1 are consistent with the values generally observed 
for AGN/GBHBs in the thermal state. For XMM3, they 
are larger and thus not constraining. In XMM2, the up- 
per limits of Ci/f are low compared with typical values 
for AGN and GBHBs, in the range 0.005-0.03. 

For the Lorentzia n model, we assum ed t he quality 
factor to be 2 as in iGoad et al.l (|200l and IHeil et all 
(2009). The upper limits of the integrated power of the 
Lorentzian R 2 are given in Table [4] for two Lorentzian 
centroid frequencies: 10~ 3 and 1 Hz. They roughly in- 
crease with the centroid frequency assumed. Compared 
with the ty pical values of 0.01 for R 2 seen in the GBHB 
hard state (Ivan der Klisl 12006). the upper limits we ob- 
tained for HLX-1 are high for all three observations. Due 



to poor statistics, we thus cannot place any firm con- 
straints on the spectral state using the timing analysis. 

5.3. Combined diagrams 

We show the best fit models for all the fitted spec- 
tra in Figure [9j which illustrates the evolution of the 
soft disk emission and the hard tail during spectral state 
transitions. For each Swift, XMM-Newton and Chandra 
spectra, we estimated the unabsorbed luminosity in the 
range 0.2-10 keV for the best fit, and calculated the flux 
ratio between the bands 0.3- 1 and 1-10 k e V. Th ose en- 
ergy ranges were also used bv lGodet et all (|2009l) . so the 
HID shown in Figure [10] is an updated version of their 
Figure 2. 

We plotted in this diagram values coming from dif- 
ferent instruments with independent calibrations. There 
could be systematic errors in the different calibrations, 
leading to a discrepancy of up to 20% in the normal- 
ization of the flux as it co uld be observed during simul- 
taneous observations (e.g. iTsujimoto et al.l l201lt ). We 
indeed observed a difference between the Swift XRT S7 
and Chandra ACIS-S spectra of HLX-1 (Figure [9]) taken 
over a similar period. However, for these two spectra, 
hardness ratio and luminosity estimates are consistent 
within this range of errors. Consistent trends are seen in 
both the Swift and XMM-Newton data, we are thus con- 
fident that the spectral variations detected in this work 
are real and significant. 

We also computed variability rms values and limits in 
the 0.0001-0.1 Hz range. They are reported in Table H 
and plotted on the HRD in Figure [TO] 

6. DISCUSSION 

6.1. Spectral states and transitions 

The high luminosity state of HLX-1 is associated with 
a soft emission with a dominant thermal disk emission, 
while the low luminosity state is found to be significantly 
harder and power-law-li ke. This is reminiscent of the be - 
havior of GBHBs (e.g. iRemillard fc McClintockl 120061 ). 
and the HID and HRD shown in Figure [TU] a re con- 
sistent with similar diagrams for GBHBs (e.g. IBellonil 
I2010D . HLX-1 therefore clearly showed transitions from 
one state to another, confirming the first evidence found 
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Table 4 

Results of the PDS Fitting 



Obs 


c P 


xi M 


Ci/f (io- 


4 , BPL) 


R 2 (10- 


' 2 , BLN) 


rms(%) 








10" 3 Hz 


1 Hz 


lO- 3 Hz 


1 Hz 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


XMM1 


46.11±0.22 


1.13(83) 


<98.8 


<58.3 


<4.2 


<326.3 


<32.3 


XMM2 


16.02±0.01 


1.28(120) 


<2.7 


<2.5 


<0.49 


<58.8 


« 4+9-1 


XMM3 


898.11±0.55 


0.79(127) 


<646.8 


<579.9 


<79.2 


<1293.7 


<96.7 



Note. — Columns: (1) Observations; (2) Best-fitting constant Poisson level in 
[rms/mean] 2 Hz -1 ; (3) Reduced \ 2 and degrees of freedom of the fits with a constant 
Poisson level; (4) Best-fitting C\ if assuming a break frequency of 10 — 3 Hz; (5) Best-fitting 
C\i r assuming a break frequency of 1 Hz; (6) Best-fitting R 2 assuming a Lorentzian ccn- 
troid frequency of 10 -3 Hz; (7) Best-fitting R? assuming a Lorentzian centroid frequency 
of 1 Hz; (8) Fractional rms variability within 0.0001-0.1 Hz, assuming the Poisson level 
to be the average PDS above 1 Hz. All limits are at a 90% confidence level, except for 
the fractional rms (column 8), which is la. 
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Figure 9. Models of X-ray spectra of HLX-1 in different states as 
obtained from the best fit to the data (see Table [2}. For XMM3, 
the abs(pow)+abs(mekal) model is shown as a dashed line, and the 
model without the mekal component as a thick line. For Chandra, 
the effect of pile-up has been corrected. 



bv lGodet et al.l (|2009f ). However, fewer luminosity states 
have been observed and therefore the diagrams are not 
equally sampled than for GBHBs. 

The high luminosity state of HLX-1 is well described 
by an optically thick disk model with temperatures of 
0.18 to 0.26 keV (Table [2]). Moreover, the low variabil- 
ity level in the XMM2 data (~10%) can reasonably be 
ascribed to the dominance of the thermal disk compo- 
nent. All this is consistent with HLX-1 being in a ther- 
mal state. The temperature of the thermal component is 
4 to 5 times lower than for typi cal GBHBs with stellar- 
mass black holes (~ 1 keV, e.g. lRemillard fc McClintockl 
2006). If this soft component we measured is due to 
emission from the inner region of an accretion disk, and 
the disk extends close to the ISCO, then the temperature 
is related to the mass of the black hole by the relation 

Tin oc MgH /4 (|Makishima et al.ll2000D . It is thus possible 
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Figure 10. Hardness-intensity (top) and hardness-rms diagrams 
(bottom) for HLX-1. This plot includes data from Swift, XMM- 
Newton and Chandra. The hardness ratios and luminosities are 
model dependent, and were obtained from the best fit of the spec- 
tra or using a conversion factor for the Swift XRT. XMM3' cor- 
responds to the model with the mekal component removed (see 
Section [4j . We show part of the HID with linear axes as an inset 
to better follow the evolution of the source in the high state. 

that the lower temperature is due to the presence of a 
black hole with a higher mass — by a factor of ~ 200 to 
500 — than a typical stellar-mass black hole (~ 10 M@), 
leading to a possible mass of few 10 3 M P >. This argu- 
ment has already been used bv lMiller et al.l (|2004l ) for 6 
bright ULXs (> 10 40 erg s _1 ) to strengthen their classi- 
fication as IMBH candidates even if the thermal compo- 
nent was not dominant. In our observations of HLX-1 in 
the thermal state, this component is clearly dominant. 
A complementary stu dy of HLX-1 in the thermal state 
bv lDavis etafl (201 ll) using the relativistic disk model 
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Figure 11. Disk bolometric unabsorbed luminosity as a function 
of disk temperature kT ln (see values in Tabled- The best fit with 
a fixed Ldisk oc TJ* is indicated with a dashed line. The XMM3 
point should be taken with caution as the disk component cannot 
be confirmed nor ruled out in the X-ray spectrum. 



BHSPEC provided robust lower and upper limits with 
3000 < M < 3 x 10 5 M©, firmly placing HLX-1 in the 
IMBH regime. Godet et al. (2011, in preparation) tested 
additional models to physically constrain the nature of 
HLX-1 and they also found mass estimates in the IMBH 
range. 

We report in Figure [11] the temperature and bolo- 
metric luminosity of the disk in the thermal state 
from the S3, S7, Chandra, XMM2 observations (see 
also Godet et al. 2011, in preparation). All those 
measurements in the thermal state are consistent with 
the correlation between luminosity and temperature ex- 
pected from a geome trically thin, optically thic k accre- 
tion disk (L oc T 4 , iShakura fe Sunvaevlll973T) . This 
supports the idea of HLX-1 being in a generic ther- 
mal state as this relatio n is also observed fo r GB- 
HBs in this state (e.g. IKubota fc Makishimal I2004t 
IRemillard fc McClintockl I2006D . It is however d ifferent 
to the ULX branch reported e.g. by ISorial (|2007T ) . where 
an anti-correlation between disk luminosity and temper- 
ature is seen for two ULXs. 

During the XMM3 observation at low luminosities, the 
spectrum of HLX-1 is mainly described by a power law 
model with a photon index that falls in the range 1.4 
to 2.1 (Table [5]), indicating a hard state as observed for 
GBHBs (we presented preliminary results in lFarrell et al.l 
1201 1|) . We could not constrain the rms variability of the 
source during the XMM3 observation due to low statis- 
tics. Therefore, a high rms value of ~30%, as observed 
for GBHBs in the hard state fe.g. IBellonill2010T ). cannot 
be ruled out and would be consistent with our data. As- 
suming that this hard state is similar to the one observed 
for GBHBs, one w ould expect enhanced radio emission 
(|Fender et al.l l2004). Such a detection would further con- 
firm this state as the canonical hard state of GBHBs and 



would allow a mass e stimate using the fundamental plane 
(jMerloni et al.|[200l . 

In the XMM3 spectrum, dominated by a power law 
component, we found marginal evidence for a soft disk 
component with temperature kT ln = 0.07 ± 0.04 keV, 
with a significance of 93%. Such a disk component would 
be consistent with the L oc T 4 relation reported in Fig- 
ure [Til as the expected temperature for a disk unab- 
sorbed luminosity of 4 x 10 40 erg s~ 4 is 0.09 keV. A 
disk could thus be either extending down to the ISCO, 
or maybe truncated as the best fit temperature value is 
lower than the temperature expected from the L oc T 4 
relation. This would again suggest a resemblance with 
GBHBs. 

During the XMM1 observation of HLX-1, the spec- 
trum is found to be well fitted with a steep power law 
of photon index ~3.5. Such a spectral state is also 
observed for some GBHBs (st eep power law state, e.g. 
IRemillard fe McClintockl [2006D . supporting the similari- 
ties between HLX-1 and GBHBs. 

The HLX-1 long term variability seen in Figure [TJ 
presents similarities with some lightcurves of G BHB (e.g. 
IGierlinski fc Newton! 120061: IDone et al.l 120071 ) and more 
generall y of soft X-ray tran sients when they show out- 
bursts (|King fe Ritterl [l998h . We observed a fast rise 
(^10 days or less), followed by an exponential decay (100 
to 200 days) of the X-ray flux from HLX-1 (FRED-like 
outburst). This could be explained by theoretical mod- 
els of th e hydrogen ionizatio n instability controlling the 
fast rise (King & RitteH [l998l ). The cause of this variabil- 
ity ha s been investigated in more details bv lLasota et al.l 
(|2011l ). They concluded that HLX-1 is unlikely to be ex- 
plained by a model in which outbursts are trigg ered by 
therm al- viscous instabilities in an accretion disc (jLasotal 
120011) . and they argue that a more likely explanation is 
a modulated mass-transfer due to tidal stripping of a 
star in an eccentric orbit around the IMBH. This would 
still explain a change in the accretion rate, leading to 
the state transitions we observe, and might also explain 
the rarity of a source like HLX-1, which shows unique 
properties for a ULX. 

6.2. Comparison with other ULXs 

It was sugg ested that most ULXs (few 10 39 to 
10 41 erg s _1 , with power law spectra) are stellar-mass 
systems accreting at Eddington ratios of the order o f 
10-30, with mild beaming factors b > 0.1 (King 2009). 
If a supercritical accretion disk — as proposed in the 
Galactic system SS433 for example — is seen face-on, 
the expected lumino sity could indeed reach 10 41 erg s _1 
(|Fabrika et al.ll2007t) . In the case of HLX-1 in the highest 
state, the Eddington ratio w ould reach 1 7 and the beam- 
ing fa ctor b ~ 2.5 x 10" 3 (jKined I201lh . IFreeland et~al~l 
(2006) calculated the broad-band radio-X-ray spectra 
predicted by micro-blazar and micro-quasar models for 
ULXs. They argued that a disk and a jet could be present 
in the system in a high/hard state close to the Eddington 
luminosity with a high accretion rate and gravitational 
energy released all the way down to the ISCO, leading 
to a hard spectrum (r=l. 4-2.1) at high observed lumi- 
nosities (> 4 x 10 39 erg s _1 ). They also proposed the 
existence of milli-blazars (IMBH and beamed emission) 
with luminosities > 10 erg s . The spectrum of HLX- 
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1 in the low luminosity state might be explained by such 
a model. However, the transition to a thermal state at 
even higher luminosities is not consistent with this pic- 
ture. The increase in luminosity and softening of the 
source (decrease of the peak energy as can be seen in 
Figure [9]) are not consistent with relativistic be aming 
varia tions which would have an opposite effect (e.g. lKina 
2009). All this indicates that the spectral variability we 
o bserve for HLX-1 i s not due to beaming variations. 

iGladstone et al.l (|2009f ) showed that simple spectral 
models commonly used for the analysis and interpreta- 
tion of some ULXs (power law continuum and multi- 
color disk black body models) are inadequate for a small 
sample (12 ULXs) of nearby, low luminosity ULXs that 
have high quality spectra. Two near ubiquitous features 
are found in the spectra: a soft excess and a rollover in 
the spectrum at energies above 3 keV. They suggested 
the existence of a new ultra-luminous state with super- 
Eddington accretion flows. This would favor the presence 
of stellar-mass black holes rather than IMBHs for some 
ULXs. In the spectra of HLX-1, which reaches lumi- 
nosities two orders of magnitude higher than the average 
ULX, we found no evidence of such a break above 3 keV. 
However, the number of counts, even in XMM2, is prob- 
ably insufficient to detect this feature ([Gladston e et al.l 
|2009jclaim the need of 10 000 counts). Godet et al. (2011, 
in preparation) discuss in more details the comparison 
with ULXs claimed to be in a ultra-luminous state, but 
the variety of X-ray spectra we observe for HLX-1 in this 
work, all in the ultra-luminous range, is clearly unique 
for a ULX. In particular, the thermal state observed at 
high luminosities (therefore probably at the highest ac- 
cretion rate) for HLX-1 is clearly distinct to the proposed 
ultra-luminous state spectrum for some ULXs. 

Finally, the X-ray observations reported in this work 
indicate that HLX-1 more resembles GBHBs scaled to 
higher luminosities (three orders of magnitude) than any 
other category of ULXs. In this picture, the bolometric 
luminosity (hence the X-ray luminosity) of the source in 
the thermal state should be lower than the Eddington lu- 
minosity for the system, giving a limit on the mass of the 
black hole of >9 000 M Q . The unabsorbed luminosity in 
the hard state is 2.2% of the highest luminosity recorded 
in this work for which we assume a sub-Eddington lumi- 
nosity (see Table[2]). This level of luminosity is consistent 
with the fraction of Eddington lu minosity observed for 
GBHBs in the hard state (<2% e.g. iMaccarone fc Coppl 
120031 : [Gierlihski fc Newtonll2006t ). further supporting the 
analogy of HLX-1 with GBHBs. It is generally observed 
for GBHBs that in the hard state and the steep power law 
state, the X-ray lu minosity is 1% and 100% of the Ed- 
dington luminosity ([Remillard fc McClintockl 120061 ) . re- 
spectively, leading to complementary mass estimates of 
16 000 and 20 000 M . 

7. CONCLUSION 

We observed the bright ULX ESO 243-49 HLX-1 with 
different X-ray observatories at luminosities ranging from 
2.1 x 10 40 to 1.25 x 10 42 erg s -1 . The spectral analy- 
sis revealed that HLX-1 showed three states: a thermal 
state at high luminosity, a hard state at low luminosi- 
ties, and a steep power law state. These states are typ- 
ically observed for GBHBs, which have luminosities at 
least three orders of magnitude lower than HLX-1. The 



clear spectral changes we observed are inconsistent with 
models of milli-blazars where the emission is strongly 
beamed toward us, and the spectra of HLX-1 does not 
resemble spectra of ULXs proposed to be in an ultra- 
luminous state with super-Eddington accretion. Instead, 
the source is much better explained as being homologous 
to GBHBs, which have non-beamed and sub-Eddington 
emission. The relatively low temperature of HLX-1 in 
the thermal state suggests the presence of an IMBH of 
few 10 3 Mq. Finally, in this picture, the source lumi- 
nosity should be lower than the Eddington luminosity, 
which converts to a lower limit on the mass of the black 
hole of >9000 M , placing HLX-1 in the IMBH range. 
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